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Abstract 

We investigate an S = 1/2 two-leg spin ladder with a cyclic four-spin exchange 
interaction whose interaction constant is denoted by J4, by using the density ma- 
trix renormalization group method. The interchain and the intrachain interaction 
constant are denoted by Jmng and Jieg, respectively and assumed to be antiferro- 
magnetic. It turns out that a spin gap between the singlet (S'fo^ = 0) and the triplet 
i^tot = 1) states vanishes at J^/Jicg — 0.3 for Jrung = >^ieg- This result is in contrast 
with the fact that the S = 1/2 antiferromagnetic Heisenberg ladder, that is the case 
of Jrung 7^ Jiog) J4 = 0, has a spin gap for all nonzero value of interchain interaction 
Jrung > 0. We find a larger value of the correlation length for the spin-pair correlation 
function than a linear size L of the system at J^/Ji^g = 0.3 and Jmng = Jieg^ the 
correlation length ^ is about 204 times of the lattice constant for L = 84 for these 
values of interactions. We also find that the string correlation function decays rather 
algebraically than exponentially at J^/Ji^g = 0.3 and Jmng = Jieg- These results sug- 
gest that there is a quantum phase transition at J^/Jicg — 0.3 for Jrung = Jieg- We 
estimate a phase boundary where the spin gap vanishes in a J4/ Jieg — Jrung/ Jieg plane 
and obtain a consistent result with that by a perturbation theory for Jmng/ Jieg > 1- 
PACS. 05.10CC, 73.43.Nq, 75.10.Jm 
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1. Introduction 



In recent years, S = 1/2 two- leg spin ladders have attracted considerable interests from both 
experimental and theoretical points of view. The two-leg spin ladders are ideal models for 
quasi-one-dimensional materials such as SrCu203llll, for which a spin gap has been observed. 
Dagotto et al^ and Rice et al% suggested moreover that hole doping to the two-leg spin 
ladders brings the superconductivity. With respect to theoretical interest, Barnes et a/.i 
suggested that the antiferromagnetic Heisenberg ladder with S = 1/2 spins (AFHL) as 
shown in Fig.|I| has a spin gap for all nonzero values of interchain interactions Jrung > 0. 
The spin-pair correlation function of the AFHL decays exponentially, in contradiction to 
that of the S = 1/2 antiferromagnetic Heisenberg chain, which has a gapless spectrum and 
a power-law decay of the spin-pair correlation function. It is believed that the resonating 
valence bond (RVB) picture is valid for the ground state of the AFHlBI210. 

We have another low-dimensional gapped system, that is the S = 1 Heisenberg (Haldane) 
chain!. It is known as Haldane's conjecture that the antiferromagnetic Heisenberg chains 
with integral spins are gaped, while the chains with half-integral spins are gaples J. The 
valence bond solid (VBS) state of Afileck-Kennedy-Lieb-Tasaki (AKLT) modeli is believed 
to be an ideal example of the Haldane state of the S = 1 system. This exactly soluble AKLT 
model has biquadratic term —\{Si ■ SjY in addition to the usual Heisenberg Hamiltonian. 
Whiteini presented a numerical evidence for equivalence of the VBS (Haldane) state and the 
RVB state. This implies that the spin gap in both the Haldane chain and the AFHL is due 
to the same mechanism. Kolezhuk and MikeskaS also argued the equivalence of the S = 1 
Haldane chain and an 5 = 1/2 ladder which includes biquadratic terms in addition to the 
usual bilinear terms. They used a method of the matrix-product (MP) wave functions in 
order to discuss exact ground states of these systems and obtained phase transition points 
where the spin gaps remain finit(ili. 

For the S = 1/2 two-leg spin ladder, it has been suggested that frustration brings a phase 
boundary where the spin gap vanishes0'lli. The frustration is due to bilinear terms which 
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consist of the next-nearest-neighbor spins in the ladder. It is known that a combination 
of bihnear and biquadratic terms can provide a cychc four-spin exchange interactionlll'lil, 
whose interaction constant is denoted by J4. Brehmer et a/.0 pointed out that a moderate 



value of J4 for Jrung = Jieg is consistent with experimental observationgi^rEll. Furthermore 
they found that the cyclic four-spin exchange interaction reduces an amount of the spin gap 
substantially. 

The purpose of the present study is to clarify an effect of the cyclic four-spin exchange 
interaction, which consists of bilinear and biquadratic terms of spin operators and gives a 
frustration on the ladder, on the quantum phase transition. We carry out the density matrix 
renormalization group (DMRG) method for the S = 1/2 two-leg spin ladder with the cyclic 
four-spin exchange interaction and investigate the spin gap between the singlet and the 
triplet states. We find that the spin gap vanishes at J^/ Jieg — 0.3 for Jrung = Jieg, although 
we have finite spin gaps for J4/ Jicg < 0.3. Furthermore at J^/ J\cg = 0.3 for Jmng = J\eg we 
observe a larger value of the correlation length for the spin-pair correlation function than the 
system size L by assuming an exponential decay of the spin-pair correlation function. We 
also find that the string correlation function decays rather algebraically than exponentially 
at this point. Those results suggest that there is a quantum phase transition at this point. A 
phase boundary in the J^/Jieg — Jrung/ Jcg plane obtained in the present study is consistent 
with that obtained by a perturbation theory for Jrung/ Jicg > 

In section 0, we express the S = 1/2 two-leg spin ladder model with the cyclic four-spin 
exchange interaction. In section |^, we present our results by the DMRG method: we discuss 
the spin gap in subsection |3 A| , the spin-pair correlation function in subsection |3 B| , the string 
correlation function in subsection |3 C| and the phase boundary in J4/ J\eg — Jrung/ Jieg plane 
in subsection |3 D|. In section H, we give conclusions of the present study. 
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2. Two-leg spin ladder of 5 = 1/2 with a cyclic four-spin exchange interaction 



The antiferromagnetic Heisenberg ladder of S" = 1/2 (AFHL) is described by the following 
Hamiltonian: 

i i 

where Si^i and Si^2 express the Pauli spin operators on chain 1 and 2, respectively. Those 
operators are shown by circles in Fig.|I|. The intrachain and the interchain interaction con- 
stant are denoted by Jieg and Jrung, respectively, and are assumed to be antiferromagnetic: 
Jieg > 0, Jrung > 0. They are shown by solid and broken lines in Fig.|I], respectively. We 
define a parameter a by the ratio of these interactions as follows: 

CL Jrung/ Jeg- (2) 

The Hamiltonian with the cyclic four-spin exchange interaction is described as follows: 

n = T^AFHL + Ja Y.iP^,^ + (3) 

i 

where P4^i means a permutation operator of four spins. This operator expresses a cyclic per- 
mutation of four spins clockwise: (1, i) ^ + 1) (2, i + 1) ^ (2, i) (1, i) (see Fig.||). 
The inverse operator P^l expresses the permutation of those four spins counterclockwise. 
The sum of these permutation operators is expressed by a sum of spin operators as follows: 

-P4,i+-P4,/ = 

4(>S'i,i ■ Si^i^i){Si^2 ■ 'S'j+1,2) + 4(<S'j^i ■ S'j^2)('S'i+i,i ■ Si+1^2) ~ '^{^1,1 ■ 'S'j+i,2)('S'j+i,i ■ Si^2) 

+ ■ + Si^2 ■ + Si^i ■ Si^2 + -^4+1,1 ■ 'S'j+1,2) 

+1/16. (4) 

In this way, the sum of the permutation operators is expressed by a combination of bilinear 
terms and biquadratic terms of intrachain spins and interchain spins. 



3. Results by the DMRG method 



In the present study, we use the infinite system algorithm of the DMRG method. All of the 
data calculated in the present study are obtained by using the open boundary conditions. 
By dividing the present section into 4 subsections, we present our results as for the spin gap, 
the correlation length of spin-pair correlation function, the string correlation function and 
the phase diagram in turn. 



A. The spin gap 

We define the spin gap as follows: 

A(L) =EoiL,l)-Eo{L,0), (5) 

where Eq{L, 0) and Eq{L, 1) are the lowest energies for which the total value of 2;-component 
of the spin operator, namely S^^^, is or 1, respectively, for the ladder of length L; there are 
2L sites in the ladder. We show size dependences of A(L)/Jieg for J^/ Jicg = and a = 0.2 
as a function of 1/L in Fig.^. The values of m in Fig.|^ show the number of eigenstates of 
the density matrix which are kept in the DMRG method. Difference between the values of 
the spin gap calculated for m = 64 and 128 is within the size of symbols plotted in Fig.|^. 
We estimate the value of the spin gap in the thermodynamic limit L ^ oo using the values 
for m = 128. In this extrapolation, the values for 10 < L < 84 are used in order to avoid 
the effect of small system size; we obtain the spin gap A(oo)/Jieg — 0.102. 

In Fig.H, we show the extrapolated values of the spin gap A(oo)/Jicg for J^/Jieg = as a 
function of a. The case of a = 1 and J4/ Jicg = corresponds to the AFHL with Jrung = <^ieg- 
On the other hand, the case of a = corresponds to two independent chains since there is 
no interchain interaction. We do not find a gap less region for a > 0. The present result by 
the DMRG method is consistent with the assertion given by Barnes et a/.i. In the standard 
analysis of experimental data for neutron scattering, nuclear magnetic resonance (NMR) 
and nuclear quadrupole resonance (NQR)liHil, a ~ 0.5 is expected by assuming J4 = 0. We 



obtain the spin gap A(oo)/Jicg — 0.21 for J^/ J\cg = and a = 0.5. On the other hand, 
Ja/ J\cg is estimated to be 0.07 from comparison between experimental data and numerical 
results given in Ref.|I^ In Fig.^, the spin gap for Ji/Jicg = 0.07 and a = 1 is shown as a 
function of 1/L. The behavior in Fig.^ is a typical example of size dependence of the spin 
gap in the DMRG method, which occurs for the system with a finite value of the spin gap. 
We estimate the spin gap for J^/Jicg = 0.07 and a = 1 to be 0.218 in the thermodynamic 
limit. We notice that the both cases, a = 1, Jij J\cg = 0.07 and a = 0.5, J^/ Jieg = give 
similar values for the spin gap. Hence there is a possibility that the assumption J4 = is 
not valid for the analysis of results for NMR or NQR. 

In our calculations, the spin gap decreases as the value of J4 increases. We show the 
result for J4/ Jieg = 0.3 and a = 1 in Fig.^. We find that the values of the spin gap for L = 4n 
and for L = 4n + 2 (n = 1, 2, ■ ■ ■) show different behavior as for size dependence, especially 
for smaller system sizes. Therefore, the data for L < 20 are discarded for an extrapolation 
to L ^ 00. Hence we obtain A(cxD)/Jicg = 0.00 ± 0.003 for J^/ J\cg = 0.3 and a = 1. 

We show the value of spin gap as a function of Ji/ J\cg for a = 1 in Fig.|^ As the value 
of JaI J\eg increases, the value of spin gap decreases gradually. We notice that the spin gap 
has very small values in the region J4/ Jieg ^ 0.2 and becomes zero near J4/ Jieg ~ 0.3. 



B. The Correlation length of the spin-pair correlation function 

We estimate the correlation length of the spin-pair correlation function {Sf ^S^i). We choose 
sites i and j such that they are located symmetrically with respect to the center of the ladder 
as shown in Fig.||. We have obtained the same results for {Sl2Sj 2) as those for (5'fiS'|j), 
and hence we show the results only for In FigJ^, the spin-pair correlation function 

for J4/ Jieg = 0.25 and a = 1 is shown as a function of distance between spins as a typical 
example. For the gapfull case, we assume the asymptotic behavior of the spin-pair correlation 
function as follows: 

(Sl.Sl,) oc e-^. (6) 
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On the other hand, we expect the power-law decay of the spin-pair correlation function in 
the gapless case as follows: 

(^•i^^i) oc \^-Jr, (7) 

where r) means the critical exponent for the correlation function. Inset of Fig.|^ is a semilog- 
arithmic plot of the absolute value of the spin-pair correlation function. A deviation from 
linear shape around |« — j | ~ 80 is due to the effect of the open boundary. We have used 
a range 30 < L < 70 for estimation of the correlation length ^. The correlation length ^ is 
estimated to be 58.6 for J4/ Jieg = 0.25 and a = 1 by use of the least squares method for the 
semilogarithmic plot. This value of the correlation length is much larger than ^ ^ 3.2 for 
the AFHL, namely Ji/J\eg = and a = 1. In the fitting by the least squares method, we 
obtain the error sum of squares (ESS) is 0.00029 for J4,/ J\cg = 0.25 and a = 1 by assuming 
the form (P). On the other hand, we obtain the ESS is 0.0034 by assuming the form with 
7] = 0.814 as the best fitting. These results suggest that the spin-pair correlation function 
decays exponentially. 



Next, we show the spin-pair correlation function for J4/ J\eg = 0.3 and a = 1 in Fig.|TO. 
We have obtained a quite large value for the correlation length, that is ^ = 203. 6, by assuming 
the exponential decay (|D where the ESS is 0.0058. By using the form (|^), we obtain 0.0045 
for the ESS with 77 = 0.242. Then, there is a possibility that the spin-pair correlation 
function decays in the power-law in this case. These results are consistent with the result 
that there is no spin gap for this value of J4/ Jieg. We notice that the value of 77 ~ | for this 
quantum phase transition is different from ?7 = lforS' = l/2 spin chainiii. 

C. The string correlation function 



The string correlation function is defined as follows: 

9{\^-j\) = (snU e-'^S]) , (8) 



where 



We notice that S'f in eq.(H) consists of two S = 1/2 spins as defined by eq.(|). These two 



spins are also illustrated in Fig.|Tl]. This choice of two S = 1/2 spins was also adopted for the 
antiferromagnetic two-leg ladder without J4 by White0, because this pair of two 5 = 1/2 
spins is expected to become effectively a single S = 1 spin in the antiferromagnetic two-leg 
ladder. We choose sites i and j such that they are located symmetrically with respect to 
the center of the ladder as shown in Fig.H as well as for the spin-pair correlation function. 

In Fig.|l2|(a), we show values of the string correlation function for some values of J^/Jieg 
and a = 1 as a function of distance |« — j| between effective spins, S'f and Sj. For J4 = 0, 
namely for the AFHL, the value of \g{\i — j\)\ takes about 0.3801 in range 5 ^ \i — j\ ^ 25, 
where there is no effect of a short distance between and and the open boundary. This 



value of Is'di— is close to the value of 1(7(0)0)1 = 0.38010765 which was obtained in Ref.|Tl 
The value of \g{\i — decreases as the value of J4 increases. We have effect of a short 
distance due to small values of |2 — j| and that of open boundary for large values of — It 
should be noticed that these effects increase as the value of J4 increases. Decreasing behavior 
of Ififd^ — j|)| is shown in Fig.[l^(b) as a function of J4/ Jieg by using values of | (7(1 7) | in order 
to avoid these two effects. This decreasing behavior suggests that the string correlation 
function can vanish for J4,/ J\cg ^ 0.3. 

We investigate an asymptotic behavior of a decay of the string correlation function at 
a = 1 and J^/Jieg = 0.3. We carry out the infinite algorithm of the DMRG method until 
a system size increases up to L = 82. We use a range of system size 7 < L < 25 in order 
to avoid the short distance effect and the open boundary effect. In Fig.|13|(a), values of 
log{|5'(|2 — J 1)1} are shown as a function of \i — j\. If the string correlation function decays 
exponentially as the increasing distance, we should have a straight line in this figure. In 
Fig.|T3|(b), values of log{|(7(|z — j|)|} are shown as a function of log(|i — j|). If the string 
correlation function decays algebraically as the increasing distance, we should have a straight 
line in this figure. Comparing Fig.|T3|(a) and Fig.^(b), we find that a line in Fig.|T3|(b) 
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approximates to a straight line better than that in Fig.|T3|(a). Hence, the string correlation 
function decays algebraically rather than exponentially at a = 1, J^/Jicg = 0.3. This result 
suggests that a = 1, J^/Jicg = 0.3 is a critical point of a quantum phase transition. 



D. Phase diagram 

We estimate a phase boundary by searching for values of J4/ Jicg at which the spin gap A(oo) 
vanishes for fixed values of a. The obtained phase boundary is shown in Fig.0 in J4/ J\eg — a 
plane, where a = Jmngl J\cg- In this J^/Jieg — d plane, the point (0,1) corresponds to the 
AFHL which has a spin gap and a point (0,0) to the antiferromagnetic chain which does not 
have a spin gap. As a decreases, the critical value of J4 decreases. In the present study, we 
are not able to conclude whether we have a finite range of J4 with a finite spin gap at a = 
or not. 



A dotted line in Fig.|T^ shows a phase boundary estimated from the spectrum of the lowest 
triplet excitation by Brehmer et a/.0 who obtained it by a perturbation theory assuming 
</icg/-^mng ^ 1 and Ji/ Jrung ^1. In the limit Jieg = and J4 = 0, the singlet dimers are 
located every rung of the ladder. Their result was obtained by including terms of third order 
in Jieg and J4. Although we can not compare the results by the perturbation theory with our 
results by the DMRG method near the region where Jmng/ Jieg — 1 and 0.25 ^ J4/ Jieg ^ 0.35, 
both results indicate that the gapfull region broadens as the value of Jrung/ Jeg increases. 

In a region where the value of J^/Jieg is lager than that of the phase boundary, the 
DMRG method becomes unstable. Hence it remains an open question whether we have a 
finite spin gap or not in that region. 



4. Conclusion 

We have investigated a quantum phase transition for the antiferromagnetic Heisenberg ladder 
with the cyclic four-spin exchange interaction by using the DMRG method; its interaction 
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constant is denoted by J4. The infinite algorithm with open boundary conditions has been 
used in order to calculate the spin gap, the spin-pair correlation function and the string 
correlation function. For J4 = 0, we have not found the gap less region by varying the 
value of a = Jmng/ Jicg- This result is consistent with the assertion given by Barnes et a/.i; 
the AFHL has a spin gap for all nonzero interchain interaction. On the other hand in the 
case of a = 1 and J4 > 0, we have found that the spin gap vanishes at Ji/Jieg — 0.3. At 
this point, a larger value of correlation length than the system size is found by assuming 
exponential decay of the spin-pair correlation function. We have obtained a better fitting 
by assuming a power-law decay of the spin-pair correlation function. We have found that 
the string correlation function decays algebraically rather than exponentially at this point. 
These results suggest that there is the spin gap for J^/Jieg ^ 0.3 and we have a quantum 
phase transition at J^/Jieg — 0.3 when Jrung = ■^'icg- This is a contrast to the result that the 
spin gap remains finite even at phase transition points for the system with a combination of 
biquadratic terms, which are not related to the cyclic four spin exchange interactional. 

We have estimated the phase boundary for the J^/Jieg — a plane and found that it is 
consistent with the result by the perturbation theory for Jrung/ Jicg > 10. Although we have 
obtained 77 ~ | at a = 1 and Ji/ J\cg = 0.3, estimation of critical exponents remains as a 
future work in order to argue the type of the critical behavior. 

The authors wish to thank T. Sakai and H. Yokoyama for fruitful discussions. This work 
was supported by the Grant-in- Aid for Science Research from the Ministry of Education, 
Science and Culture(11780183). 
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FIGURES 
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Fig. 1. Antiferromagnetic Heisenberg ladder with S = 1/2 spin (AFHL). The Pauli spin operators 
are denoted by circles. Antiferromagnetic interactions Jicg along each chain and Jrung between the 
chains are denoted by solid and broken lines, respectively. 
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Fig. 2. Antiferromagnetic Heisenberg ladder with a cyclic four-spin exchange interaction J4. 
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Fig. 3. 



A(L)/Jieg as a function of 1/L for Ji/J\eg = and a = 0.2. 




Fig. 4. Extrapolated value of spin gap for J^/Jieg = as a function of a. 
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Fig. 5. Size dependence of the spin gap for Ja/ J\eg = 0.07 and a = 1. 
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Fig. 6. Size dependence of the spin gap for J^/ Jieg = 0.3 and a = 1. 
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Fig. 7. J4/Jieg dependence of the spin gap A(oo)/Jieg for a = 1. 
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Fig. 8. A choice of sites i and j. We assume that they are located symmetricany with respect 
to the center of the ladder. 
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Fig. 9. Spin-pair correlation function as a function of distance between Sf i and Sj i for 
Ji/ J\eg = 0-25 and a = 1. 
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Fig. 10. Spin-pair correlation function as a function of distance between two spins for 

Ji/ J\eg = 0-3 and a = 1. 



Fig. 11. A single effective S = 1 spin which consists of two 5 = 1/2 spins enclosed by dotted 
line for the antiferromagnetic two-leg spin ladder. 
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Fig. 12. (a) String correlation function for cases in which the value of J^/ J\f,^ is 0.0, 0.1, 0.2 
and 0.3 as a function of distance \i — j\ between Sf and S^. (b) Decreasing behavior of the string 
correlation function for |i — jl = 17 as a function of J4/Jieg. Both figures are for the case a = 1. 
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Fig. 13. (a) The value of log{|y(|i — j|)|} as a function of |i — j|. (b) The value of log{|(7(|z — j|)|} 
as a function of log(|i — 
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